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The Hawking radiation is one of the most interesting phenomena predicted by the theory of 
quantum field in curved space. The origin of Hawking radiation is closely related to the fact that a 
particle which marginally escapes from collapsing into a black hole is observed at the future infinity 
with infinitely large redshift. In other words, such a particle had a very high frequency when it was 
- near the event horizon. Motivated by the possibility that the property of Hawking radiation may be 

altered by some unknowned physics which may exist beyond some critical scale, Unruh proposed a 
model which has higher order spatial derivative terms. In his model, the effects of unknown physics 
are modeled so as to be suppressed for the waves with a wavelength much longer than the critical 
scale, fcg ^ . Surprisingly, it was shown that the thermal spectrum is recovered for such modified 
models. To introduce such higher order spatial derivative terms, the Lorentz invariance must be 
violated because one special spatial direction needs to be chosen. In previous works, the rest frame 
of freely-falling observers was employed as this special reference frame. Here we give an extension 
by allowing a more general choice of the reference frame. Developing the method taken by Corley, 
' we show that the resulting spectrum of created particles again becomes the thermal one at the 

, Hawking temperature even if the choice of the reference frame is generalized. Using the technique of 

■ the matched asymptotic expansion, we also show that the correction to the thermal radiation stays 
' of order Uq^ or smaller when the spectrum of radiated particle around its peak is concerned. 

a^ ■ 

■ I. INTRODUCTION 

cr| 

^ , The thermal radiation from a black hole was first predicted by Hawking Q|, which phenomenon became widely 
known as the Hawking radiation. This prediction is based on quantum field theory in curved space, which is thought 
^ of as an effective theory valid for low energy physics. However, when we consider the mechanism of the Hawking 
r '~j ■ radiation, crucial role is played by wave packets which left the past null infinity with very high frequency. Such 
rS \ wave packets propagate through the collapsing body just before the horizon formed, and undergo a large redshift on 
^ . the way out to the future null infinity. Here arises one question. Can it be justified to apply quantum field theory 
" " ' in curved space, an effective theory at low energy, to the phenomenon which involves the infinitely high frequency 
regime? There may exist some unknown physics which invalidates the application of the standard quantum field 
theory in curved space ||]. 

One of such possibilities is that the spacetime may reveal its discrete nature at such high frequencies. To take 
account of the effect of possible modification of theory in the high frequency regime, Unruh proposed a simple toy 
model by a sonic analog of a black hole [^ji). In Unruh model, the dispersion relation of fields at high frequencies are 
modified so as to eliminate very short wavelength modes. In some sense, this modification is arranged to reflect the 
atomic structure of fluid which propagates sound waves. Usually, the group velocity of sound waves drops to values 
much less than the low frequency value when the wavelength becomes comparable to the atomic scale. In performing 
such modifications [Q, one must assume the existence of a reference frame because the concept of high frequencies 
can never be a Lorentz invariant one. Namely, the Unruh's model manifestly breaks the Lorentz invariance. To our 
surprise, even with such a drastic change of theory, the spectrum observed at the future infinity was turned out to be 
kept unaltered 0-0]. Here, in this paper, we consider a generalization of this model. 
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The Lorentz invariance is the very basic principle for both the special relativity and the general relativity. Hence, 
there are many efforts to examine the violation of the Lorentz invariance , and new ideas to make use of high energy 
astrophysical phenomena are also proposed recently |^. However, we have not had any evidence suggesting this rather 
radical possibility yet. Therefore, one may think that it is not fruitful to study in detail such a toy model that violates 
the Lorentz invariance at the moment. But, we also have another motivation to study this toy model even if we 
could believe that the Lorentz invariance is an exact symmetry of the universe. In most of literature, the Hawking 
radiation was studied in the framework of non-interacting quantum fields in curved space. However, when we consider 
a realistic model, it will be necessary to consider fields with interaction terms ^ . The evolution of interacting fields 
in the background that is forming a black hole is a very interesting issue but to study it is very difficult. Hence, as a 
first step, it will be interesting to take partly into account the interaction between the quantum fields and the matter 
which is forming a black hole. Then, it will be natural to introduce a modified dispersion relation associated with the 
rest frame of the matter remaining around the black hole. In this sense, the Unruh's model does not require that the 
fundamental theory itself violate the Lorentz invariance. 

In order to introduce the modified dispersion relation we need to specify one special reference frame. In previous 
works the rest frame of freely-falling observers was employed as the special reference frame. In this case, the 

thermal spectrum of the Hawking radiation was shown to be reproduced. However, it is still uncertain whether the 
same thing remains true even when we adopt another reference frame as the special reference frame. In this paper, 
we give a generalization of previous works P-|^,|lO| by allowing a more general choice of the reference frame. 

In most parts of the present paper, we follow the strategy taken in the paper by Corley ||^. (See also In 
his paper, as modifications of the Unruh's original model, two types of models were investigated. One is that with 
subluminal dispersion relation and the other is that with superluminal dispersion relation. It was shown analytically 
that the thermal spectrum at the Hawking temperature is reproduced in both cases. However, in the superluminal 
case, the standard notion of the causal structure of black hole breaks down. Even if we consider the case that the 
background geometry is given by a Schwarzschild black hole, the wave packets corresponding to the radiated particles 
can be traced back to the singularity inside the horizon due to their superluminal nature. Hence, the singularity 
becomes naked, and we confront the problem of the boundary condition at the singularity. To avoid this difficulty, it 
is often required that the vacuum fluctuations be in the ground state just inside the horizon. However, it is not clear 
what is the correct boundary condition. As a topic related to the superluminal dispersion relation, it was also reported 
that the Hawking radiation is not necessarily reproduced in the models with an inner horizon pl] |. In this paper, we 
wish to focus on the subluminal case leaving such a delicate issue related to the superluminal dispersion relation for 
the future problem. Even in the restriction to the subluminal case, it will be important to study various models to 
examine the universality of the Hawking radiation. In the present paper, we finally find that the resulting spectrum of 
created particle stays thermal one at the Hawking temperature as long as we mildly change the choice of the special 
reference frame. By a systematic use of the technique of the so-called matched asymptotic expansion, we also evaluate 
how small the leading correction to the thermal radiation is. On the other hand, for some extreme modification of 
the reference frame, in which case the analytic treatment is no longer valid, the spectrum is numerically shown to 
deviate from the thermal one significantly. 

This paper is organized as follows. First we introduce a generalization of the Unruh's model in Sec. 2. In Sec. 3 we 
review what quantities need to be evaluated in computing the spectrum of particle creation in our model. In Sec. 4, we 
construct a solution of the field equation, and we evaluate the spectrum of created particles by using this solution. To 
determine the order of the leading correction to the thermal spectrum, we employ the method of asymptotic matching 
in Sec. 5. In Sec. 5, we also demonstrate some results of numerical calculation to verify our analytic results. In addition, 
we display the results for some extreme cases which are out of the range of validity of our analytic treatment. Section 
6 is devoted to conclusion. Furthermore, appendix C is added to discuss the effect of scattering due to the modified 
dispersion relation. Although we do not think that this effect is directly related to the issue of Hawking radiation, it 
can in principle change the observed spectrum of Hawking radiation drastically if it accumulates throughout the long 
way to a distant observer. In the present paper, we use units with h ^ c = G = 1. 



II. MODEL 

In this section, we explain how we generalize the model that was investigated in the earlier works [Q-^- Following 
these references, we consider a massless scalar field propagating in a 2-dimensional spacetime given by 

ds"^ -dp + {di - v{x)dif , (2.1) 

where v{x) is a function which goes to a constant at x ~* oo and satisfies v{x) > — 1 for i > 0. The equality holds 
at i = 0. Since the line element di = is null at v — —I, we find that the event horizon is located at v = — 1. 
Furthermore, by the coordinate transformation given by dt = dt' + v/{l — v^)dx, the above metric can be rewritten as 
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ds"^ = -{1 ~ i'^)di''^ + —dx'^. (2.2) 

1 — 

If we set v{x) — —y/2M/ {x + 2Af), this metric represents a 2-dimensional counterpart of a Schwarzschild spacetime 
with the event horizon at 5; = 0. In this coordinate system, the unit vector perpendicular to the t =constant 
hypersurfaces is given by Uq, := dat, and the differentiation in this direction is given by u°'da = {di + vd^). We denote 
the unit outward pointing vector normal to Uq, by s" . In order to examine the effect on the spectrum of the Hawking 
radiation due to a modification of theory in the high frequency regime, they investigated a system defined by the 
modified action of a scalar field, 

S=\l (fiV^g'^^V^rVpcj,, (2.3) 



where the differential operator V is defined by u'^Va = u"da, s^Va = F{s°'da). If we set F{z) = z, the action ( |2.3| ) 
reduces to the standard form. Since we are interested in the effect caused by the change in the high frequency regime, 
we assume that F{z) differs from z only for large z. 

In the above model, the dispersion relation for the scalar field manifestly breaks the Lorentz invariance, and there 
is a special reference frame specified by u. One can easily show that this reference frame is associated with a set 
of freely-falling observers. As noted in Introduction, it was shown that the spectrum of the Hawking radiation is 
reproduced in this model. Here we consider a further generalization of this model allowing to adopt other reference 
frames as the special reference frame. 

However, because of technical difficulties, we restrict our consideration to stationary reference frames. As we are 
working on a 2-dimensional model, the reference frame is perfectly specified by choosing one time-like unit vector, 
which we denote by u. Since 9^ in the original coordinate system (i, x) is a time-like Killing vector, the condition for 
the reference frame to be stationary becomes ^Ojit" — 0, where £,q~ is the Lee derivative in the direction of d^. This 
condition can be simply written as du^/dt = 0, where we used indices associated with ""to represent the components 
in the (f, x)— coordinate. Furthermore, the covariant components Ua{x) = g^jj{x)u^{x) is also independent oil. Thus, 
if we introduce a new time coordinate t by 

dt = u'^^ {x) {u^{x)di + ui{x)dx') —di~^{x)dx, (2-4) 

the constant hypersurfaces become manifestly perpendicular to Ua- Here we introduced 7(5;) := —Ui{x)/uq{x). 
Further, it is convenient to choose a new spatial coordinate x so that dt coincide with the Killing vector Since 

d,=di+^d,, (2.5) 
ot 

X should be chosen as a function which depends only on x. Hence, we set 

x = f CiS:')di'. (2.6) 





By using such a new coordinate {t, x) with 

C(i') = (1 - ^^7)' - 7', (2.7) 

the metric (2J) can be written in the form|], 

ds^ = n^{x) [-dt^ + {dx - v{x)dtf] , (2.8) 



where we set 



v{x) -.^ 'f + v{l - v'^). (2.10) 



By considering sonic analog of the Hawking radiation, the use of this type of conformal metric was discussed [Q. 
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Here we mention the constraint on ft. If we explicitly write down the condition that u is a time-like unit vector, i.e., 
Mqu" = —1, we find that fi^ = u? holds. Hence il^ > is guaranteed as long as u is time- like. Also, directly from the 



metric (2.8), we can easily verify fi^ > if and only if the t— constant hypersurfaces are space-like. Hence, it will be 
appropriate to assume > 0. Then, we find that fl^ has a finite minimum value when |ti| < 1. The minimum value 
is 1 — which is realized when v = 0. 

Next we write down the explicit form of u and s. By using the facts that u is perpendicular to the t— constant 
hypersurfaces and that it is an unit vector, we can show that the differentiation in the direction of u is given by 



= ^ (dt 



vdx 



(2.11) 



Similarly, we can show that the differentiation in the direction of s, which is an unit vector perpendicular to u, is 
given by 



(2.12) 



By using u and s, the metric can be represented as g°'^ = —u°'u^ + s'^s^ , and the determinant of gap becomes 

Now, we find that to generalize the choice of the special reference frame introduced to set the modified dispersion 
relation is equivalent to generalize the metric form given in (2.1) to the one given in (2.8) replacing u and s with u 
and s in the defining equations of the differential operator V. As a result, the action ( p.3| ) becomes 



2 



dtdx 



(2.13) 



If we set = 1, the models are reduced to the original one discussed in Ref. Q. 

Here, we should mention one important relation for the later use. The temperature of the Hawking spectrum is 
determined by the surface gravity n defined by k :— dv / dx\x=o- The surface gravity k can also be represented as M] 



dv 
dx 



In order to verify this relation, we evaluate dv / dx by using Eqs. (|2j) and ( |2.10| ) to obtain 

1 



dv 
dx 



(1 — v^Y ~ T 



(2.14) 



(2.15) 



From Eq. ( 2.10 ), we also find v — —1 when a; = 0. Hence, substituting w = — 1 into ( p.l5| ), we obtain Eq. ( 2.14 ). 

Then, let us derive the field equation by taking the variation of the action ( 2.13|) . Assuming that F{z) is an odd 
bmction of z, we obtain 



(2.16) 



To proceed further calculations, we need to assume a specific dispersion relation. Following Rcf. we adopt here 

1 



(2.17) 



where fco is a constant. Since the model should be arranged to differ from the ordinary one only in the high frequency 
regime, the critical wave number is supposed to be sufficiently large. With this choice of F, neglecting the terms 
that are inversely proportional to the fourth power of fco, the field equation becomes 



di 



^ d^U U — 



(2.18) 



Before closing this section, we briefly discuss the meaning of the functions and v. From ( ^.11 ), it is easy 
to understand that v is the coordinate velocity of the integration curves of u. To understand the meaning of fJ, 
we further calculate the covariant acceleration of the integral curves of u, where semicolon represents the 

covariant differentiation. After a straightforward calculation, we see that the covariant acceleration is given by jc^j;!!"^ |. 
Hence, we find that the derivative of fi"^ gives the acceleration of the reference frame. 
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III. PARTICLE CREATION RATE 



In this section, we briefly review how to evaluate the spectrum of Hawking radiation. We clarify what quantities 
need to be calculated for this purpose. (For the details, see Ref. p.) 

To evaluate the spectrum of Hawking radiation, we need to solve the field equation ( 2.1§| ) with an appropriate 
boundary condition. However, owing to th e time translation invariance with the Killing vecto r (?t, we do not have to 
solve the partial differential equation ( 2.1^ ) directly. By setting (j){t, x) = e~'^^*'ip{x), Eq. (^.li 
differential equation (ODE), 



reduces to an ordinary 



{—iuj + dxv){—iuj + vdx) — s dx^ + 



2fcn 



2kn 



d -d -d^ 



(3.1) 



We could not make use of this simplification if we relax the restriction that the reference frame should be stationary. 

Here we note that the norm of dt is given by \dt \ — \/l — v'^Vl = \/l — v"^. Therefore, the frequency for the 

static observers who stay at a constant x(or x) is related to lo by 



ui 



[stat] 



'I 



(3.2) 



Hence, as long as v^o '■— v{x oo) is not equal to zero, aj'^**"*^ differs from cj. By looking at the metric (2.2) in 
the static chart, we find that this frequency shift is merely caused by the familiar effect due to the gravitational 
redshift. Therefore, even if we consider models with Voo 0, lu might be identified with the frequency observed at the 
hypothetical infinity where the gravitational potential is set to be zero. However, the situation is more transparent 
if we can set Vao = like the 2-dimensional black hole case. In this case, we can identify uj with w*^"*"*' without 
any ambiguity. As mentioned in Ref. there is a difficulty in evaluating the spectrum of radiation in the case of 
Voc ■= — > oo) 0. In previous models, Voo — directly implies 5oo = 0. Therefore, we could not apply the result 
to the example of a 2-dimensional black hole spacetime directly]^ On this point, in our extended model, the cases 
with Woo = can be dealt with since Vac = does not mean Vao 7^ 0. 

Let us return to the problem to solve Eq. (3.1). From the above ODE, the asymptotic solution at x ^ cx) is easily 
obtained by assuming the plane wave solution like 



'^{x) cx e 



ihx 



Substituting this form into Eq. (3T), we obtain the dispersion relation 

where f2oo is the asymptotic constant value of fJ. The quantity on the left hand side 

Lo' := — vk, 



(3.3) 



(3.4) 



(3.5) 



is related to the frequency measured by the observers in the special reference frame. In fact, this frequency divided 
by n is the factor that appear s w hen we act the operator iu"'da on a wave function ^-^i'^t-kx) _ gjj^gwn in Ref. 
two of the 4 solutions of Eq. (3.4) have large absolute values, which we denote by k±, and the other two have small 
absolute values, which we denote k±s. For each pair, one is positive an d th e other is negative. The subscript ± 
represents the signature of the solution. Then, the general solution of Eq. ( ^.4| ) at a; ^ 00 is given by a superposition 
of these plane wave solutions as 



^iki {uj)x 



(3.6) 



l=±.±s 



By such a local analysis, however, the coefficients 01(10) are not determined. To determine them, we need to find 
the solution that satisfies the boundary condition corresponding to no ingoing waves plunging into the event horizon. 
This condition is slightly different from the condition that the solution of ODE (3.1) vanishes inside the horizon. The 
latter condition is stronger than the former one because the latter one also prohibits the pure outgoing wave from 



In the model proposed in Ref. [Q, the case with Voo = can be dealt with. 
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the event horizon which may exist in the present model with the modified dispersion relation. The former condition 
is the sufficient condition to determine the wave function uniquely, while the existence of the solution that satisfies 
the latter condition is not guaranteed in general. However, once we find such a solution that satisfies the latter 
stronger condition, it is the solution that satisfies the required boundary condition. In the succeeding sections, we 
solve ODE (3.1) requiring the latter condition. 

Finally, we present the formula to evaluate the expectation value of the number of emitted particles naturally 
defined at a; — *■ oo. In spite of our generalization of models, the same derivation of the formula that is given in Rcf. |^ 
is still valid. The same arguments follow without any change, but one possible subtlety exists on the point whether the 
expression of the conserved inner product is unaltered or not. Therefore, we briefiy explain this point. The defining 
expression for the conserved inner product given in Ref. |^ is 

(01,02) = i Jdx[(t)i*{dt +vdx)(j}2 - 02(9t + vdx)(j}i*], (3.7) 

where the integration is taken over a t-constant hypersurface. Here both 0i and 02 are supposed to be solutions of 
the field equation. The constancy of this inner product is related to the invariance of the function 



£(01, 02) := (dt + 9.)0i • (dt + 9^)02 - 
under the global phase transformation 



(3.8) 



(3.9) 



Taking the differentiation of C{e'^^4)i,e^^4>2) with respect to A, we have 
d£(e*^0i,e'^02) 







—i(t) + —. + -^(-102) + 



dC 



901 dt 
dC 



902 



dC C?( — 102) 
rlih^ dt 



dt 



(3.10) 



where we used the field equation in the last equality. Eq. ( |3.10 ) proves the constancy of the inner product (3/7). Of 
course, the constancy of the inner product (3.7) can also be verified by directly calculating its i-derivative using the 
field equation. 

Now that we verified that the extension of model does not alter the expression of the inner product, we just quote 
the formula from Ref. For a wave packet which is peaked around a frequency uj, the expectation value of the 
number of created particles is given by 



N{uj) 



\u {k.)vg{k.)cl{io)\ 
\uj' {k+s)vg{k+s)cX,{uj)\' 



(3.11) 



where Vg(k) :— dui(k)/dk is the group velocity measured by a static observer. 



IV. SOLVING FIELD EQUATION 

In this section, to determine the coefficients ci, we solve the field equation ( |3.1| ) by using several approximations. 
In the region close to the horizon, we use the method of Fourier transform. The solution is found to be uniquely 
determined by imposing the boundary condition discussed in the preceding section. On the other hand, in the region 
sufficiently far from the horizon, we construct four independent solutions which become e^'^ii^)^ at x ^ oo. We use 
the WKB approximation for the two short-wavelength modes and we use the simple l/feg expansion for the other two 
long-wavelength modes. Later, we find that these two different regions of validity have an overlapping interval as long 
as ko is taken to be sufficiently large. Hence, the requirement that the solutions obtained in both regions match in 
this overlapping interval determines the coefficients c^((jj). 

Basically, our computation is an extension of that given by Corley 1^. Here we take into account the generalization 
of models discussed in Sec. 2. Furthermore, to evaluate the order of the leading deviation from the thermal spectrum, 
we shall take into account some higher order terms. At the same time, we also carefully keep counting the order of 
errors contained in our estimation. For brevity, we concentrate on the most interesting case in which uj and k are 
same order. 
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A. the case close to the horizon 



Now we want to find a solution satisfying the boundary condition that the wave function rapidly decrease in the 
horizon. Therefore, we restrict our consideration to the region close to the horizon, \x\ < xi, by choosing a sufficiently 
small xi. We introduce a parameter. 



S := max Si{x), 

|a:|<2;i,i=0,l,2 



(4.1) 



with 



do := KX, di , 02 := — x. 

K ill 



(4.2) 



Since we wish to think of (5 as a small parameter for the perturbative expansion, we require S ^ I. Then, we find 
that xi must be chosen to satisfy 



xi <C min k/ki, fli/rio) • 

We assume that v{x) and l/n{x) can be expanded around the horizon like 

v{x) ^ -1 + KX+ + 0{S^), 



(4.3) 



r^o f^i r^o 



(4.4) 



Substituting these expansions into the field equation (3J), we classify the terms into five parts according to the 
number of differentiations acting on ip. Then, keeping the leading order correction terms with respect to S in each 
part, we obtain the field equation valid in the region close to the horizon as 



1 2 



+ {2Hix - t^x^ + k\x^)iI," 



+2 [— (iti-" — k) + iyK(iuj — k) + Kx)x] '0' ~ ioj{iuj — n — h\x)^Ij = 0. 



In Corley's paper [^, the terms of 0{5^) were neglected, while we keep them in the present paper. 
Here, we introduce the momentum-space representation of the wave function, il}{s), by 



(4.5) 



il){x) 



(4.6) 



Substituting this expression into Eq. (4.5), we perform the integration by parts like 



where Si and s f are the start and the end points of integration, respectively. Note that there appear surface terms like 
the last term on the right hand side in the above example. Then, we find that the field equation in the momentum 
space is given by 



L[)p{x)] :— / ds L[il}{s)] + (boundary terms) = 0, 
Jc 



(4.7) 



where 



ms)] 



+2 [— (iw — k) — {n{iuj — k) + K^)9s] stp — iuj{iuj — k + K'fds)il>. 



(4.8) 



If we are allowed to neglect the boundary terms, we can construct a solution of Eq. (4.5) by using Eq. (4.6) from a 
solution tp{s) which satisfies 
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L[^{s)] = 0, 



(4.9) 



In the following, we solve Eq. (4.9) without taking any care about whether the boundary terms can be neglected or 
not. After wc fi nd a solution, we verify that the corresponding boundary terms are sufficiently small. 
To solve Eq. (4.9), it is convenient to introduce a new variable 



W := 



dins 



Taking S as a small parameter, we expand W as 
Then, is foimd to be given by 



and W^(i) is given by 



w^(i) =7e(w(°)). 



with 



n{W) := ~do 



2ains 2 2 
I dW 1 



















2k J 



-^^(-2ah^-2^ +2^-l)(^^) 



Here we defined 



e := 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



(4.15) 



Although we later consider the situation in which e is small, e is not small at all in the region close to the horizon. 
Substituting the explicit expression for W^"^ into TZ{W), we obtain 



lOJ ~ lUJ , , , -I 



(4.16) 



Under the condition that the boundary terms in Eq. (4.7) can be neglected, we obtain the solution '4'{x) which is valid 
up to 0{5^) as 



c 



where, 



is given by 



xf{s) = xfois) + xfi{s), 



(4.17) 



(4.18) 



xfa{s) ■.= xs + J ^W^(o) -21n 
e'^{xs)^ 



iw , , 
-1 Ins, 



(4.19) 



ds , 
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^ l{sxf ~j[So + {Si - 2S2) ^- - 1 

iuj ~ / iu \ . , -I 



(4.20) 



The boundary condition of the wave function requires that it exponentially decreases inside the horizon. In order 
to construct the wave function that satisfies this boundary condition, we must choose the contour of integration 
appropriately. We propose to adopt the contour C given in Fig. 1. This contour does not go to infinity but have end 
points (denoted by op en circles in Fig. 1) at which the absolute value of s is sufficiently large but does not exceed 
the limit given in (Al). Hence, as shown in Appendix the higher order correction does not dominate along this 
cont our. Furthermore, as is also explained in Appendix ^ with this choice of end points, the boundary terms in 
Eq. (4.7) are exponentially small, and can be neglected. 




FIG. 1. The contour for integration is chosen for the correction terms not to dominate the leading terms. In the shaded 
region, the higher order terms with respect to 5 becomes larger than the leading terms. Hence, the contour is shortened so as 
not to violate the validity of approximation. 



Next we show that ip{x) given by (4T7) is actually the solution that satisfies the required boundary condition. In 
order to evaluate the integration ( 4.17 ) analytically, we use the method of steepest descents. For this method to be 
valid, l/|xso| ^ 1 is required, where sq is the value of s at the saddle point that dominates the integral. Then, this 
requirement can be rewritten as 



\X\ > X2, 



(4.21) 



where we must choose X2 to satisfy {l/\xso\)x=x2 ^ 1- We shall see later that l/|xso| is of 0(|e|). Hence, for the first 
time at this moment, we further restrict our consideration to the region in which e is also small. In order that the 
condition > X2 to be compatible with \x\ < xi, 



2/3 



<C min(l, k'^/k'I, kQi/V,o), 



(4.22) 



is required. However, this requirement to the model parameters will not reduce the generality of our analysis so much 
because we are interested in the case that the typical length-scale for the modification of dispersion relation, k^^, is 
sufficiently small. 

As in the case of S, we introduce an expansion parameter 

e:=\i{x2)l (4.23) 

and we neglect the terms that induce the relative error of O(e^) or smaller in the amplitude of the wave function. As 
for (5, we also keep the terms up to linear order in S. Here one remark is in order. We imposed a further restriction 
|a:| > X2 to evaluate the explicit form of the solution ( 4.17 ). We stress, however, that the solution ( 4.17 ) itself is valid 
throughout the regi on < xi. 

To evaluate ( 4.17 ) by using the method of steepest descents, we need to know the value of s at saddle points which 
are determined by solving 



/'(s)^/^(s) + /{(s) 
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lUJ 
K 



+ {sxy-w 



0. 



(4.24) 



We solve this equation by assuming that the solution is given by a power series expansion with respect to S as 

s± = so± + si± + S2± ■ ■ ■ . (4-25) 

For our present purpose, it is enough to find the solution in the form of series expansion with respect to e. One 
solution of xso is of 0(1), and the integration along the path through this saddle point cannot be evaluated by the 
method of steepest descents. The other two solutions are given by 



XSo± 



I + iuj / K 



± 



e'^ 3{l + iuj/Ky 



+ 0(6^), 



(4.26) 



and 



XSi = ± 



So 5i 



^ 32 



2 V' + - 



K 



3+^)^0 

K 



K 



-3 + — I + 8^2 





( 3iw^ 






F4 5 


1 h 


V « / 



+ 0(e2). 



(4.27) 



Since |smaa;/so±| ~ \j ^ I ^ ^ Ij thcsc saddle points are contained in the region in which the expansion with respect to 
(5 is valid. In the following, to keep the notational simplicity, we abbreviate the subscript ± from so± and si± unless 
it causes any ambiguity. 



Now we evaluate the integration (4.17) by using the method of steepest descents. For the contour given in Fig. 1, 
only the saddle point 5+ dominantly contributes to the integration inside the horizon. For our present purpose, the 
formula 



i\){x) 



(-a:/"(s±))i/2 



1 



5 {xr'\s^)f ^ 1 x/W(.±) ^ 



24 (x/"(s±))3 8 (a;/"(s±))2 



(4.28) 



is accurate enough to keep the correction up to O(e^). The details of calculation to evaluate Eq. ( 4.28 ) up to 0(e,(5) 
is given in Appendix In the end, we obtain 



V'(x) « \/27rK(fcor2o) 2(_2kx)" 



■ exp ( --^kQ^o{-xfl'^ + VF + 0(e^ (5^) 



(4.29) 



where 



W = -\\\—\ ( -— 5o + — ^1- -^2 
2 J V 10 10 5 

, lUJ 3\-: / iUJ 3\-: I- 

+ '4;: + 8r°+U:-8r^ + 2^^ 



I 4:1 iLU 1 / iuj 

yls ^ ~ ^ 4 v~ 

II 3 iuj 15 f iuj 
64 A~ 16 W 



^1 



11 1 iuj 1 f iuj\ 
64 + 4^+161-' '^^ 



5 iw 3 f ILU 
16 ~ ~ ~ 4 V~ 



(4.30) 



We can immediately see that amplitude of the wave function reduces exponentially as we decrease x (as we increase 
-x). 
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FIG. 2. The deformed integration contour to evaluate tp{x) outside the horizon. In a; > region, the saddle points move to 
the neighborhood of the imaginary axis. The contours Ci and C2 are chosen to pass through these saddle points and to be 
able to evaluate the integrations along them by the steepest descents. The remaining part of the integration contour which 
goes around the branch cut is called C3. 



Next we turn to evaluate tp{x) for a; > 0. We use the method of steepest descents again. In the present case, 
the location of saddle points moves to points on the imaginary axis on the complex s— plane. The leading order 
approximation is given by 



XSo± 



£2 



(4.31) 



Therefore, to evaluate (4.17) by using the method of steepest descents, we need to deform the contour of integration. 
At this point, we must take account of the existence of a branch cut emanating from s — 0, which originates from the 
log-term in the integrand. We choose this branch cut along the negative side of the real axis. Then, deforming the 
contour so as to go through these two saddle points we find that the contour is divided into three pieces as shown in 
Fig. ||. We respectively denote them by Ci, C2 and C3. Ci and C2 are the contours passing through the saddle points 
S- and s+, respectively. Both contours have a new boundary point which is chosen to satisfy \s\ < Smax- The contour 
C3 connects these two newly introduced boundary points, going around the origin in the anti-clockwise manner. 

First we evaluate the integrations along the contours Ci and C2 by using the method of steepest descents. Just 
repeating the same calculation as in the case of x < 0, these integrations are evaluated as 



V'i,2(a;) = e^^"/^"^— ^(fco»o)-'^%/2^(2Kx)'^-^ exp f T«^V2^fcof^oa;^/^ + T^i,2 + 0{e^,5'^)] , 



(4.32) 



where 



■\/2 / 1 J 1 J 2~ 
Wi_2 = T«— + —5i - -62 

e \ 10 10 5 




Sq 



iuj 3\ - 1 - 

:r ~ o Ml + o'52 

4k 8 / 2 



41 iuj If iuj 
48 ~ ^ 4 \~ 

11 3 iLU \b f IUJ 
64 ^ 1~ ^ 16 V"k 




1 iuj 


1 


/ iuj 




4 K 


^16 


V « 


)■) 


5 








16 


K 


i( 


K 1 



(4.33) 



Next we consider the integral along the C3 contour. Here, we divide xf{s) given in Eq. (4.18) into two parts as 

xfo{s) := xs + ^—1 — — ^ Ins 

(4.34) 



6 

and we expand e^-'^^^'*-' assuming that xfi{s) is small. From the validity of such an expansion, it is required that 
|a;/i(s)| ^ 1. As for the case with large |s|, the integrand becomes exponentially small when \xs\ » 1. There we do 
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not have to mind at all even if xfi{s) becomes large negative. In the restricted region satisfying |a;s| < 1, it is easy to 
see that xfi{s) <C 1 is always guaranteed. As for the case with small |s|, we do not have to consider the situation in 
which \s\ becomes extremely small because there is no requirement on the choice of contour except for being inside 
the saddle points. For example, if we choose the contour to be \s\ ^ 1, |x/i(s)| <C 1 is guaranteed. Therefore, we find 
that it is allowed to expand e^^^'^'^^ as 1 + xfi{s) + ■ ■ ■. 

After this expansion, introducing a new variable z by e~"z :— sx, the integration along C3 is written as 



■03 (a;) 



dz {~z) 



-e-'[l + xfi 



(4.35) 



where C3 is the contour in the complex z— plane corresponding to C3. Since the integrand becomes exponentially 
small at the boundaries, we are allowed to continue the contour to 00. Then, using the integral representation of a 
gamma function, the leading term corresponding to 1 in the square brackets of Eq. (4.35) is expressed as 



^^3(0;) ~ -2sinh(7rcj/K)r(-itj/K) x*'^/". 



(4.36) 



Next, we consider the remaining terms in Eq. ( 4.35| ). Let us express xfi as xfi = ^ a„(— z)", where the coefficients 
a„ are non-dimensional constants. Then, by using the integral representation of the gamma function, we can evaluate 
the contribution from each term as 



dz{-z)"'-^-"^^''t 



-2a„ sinh ( tt 



C3 



lUJ 
K 



(4.37) 



and we find that its relative order is simply determined by the order of a„. Hence, to find the expression correct up 
to 0(e^, S^), the only term that we must keep is 



xh{s)^~'^(^ + l]~6,{sxr 



Thus, we finally obtain 



\ 2k 



(4.38) 



(4.39) 



In this section, we approximately solved Eq. (4.5) with the boundary condition that the wave function decreases 
exponentially inside the horizon. We evaluated the explicit form of the approximate solution in the region xi > x > X2 

as 

i^{x)^^i{x)+il)2{x)+^z{x), (4.40) 
where each component ^i{x) is given by Eq. (0|) or Eq. (4.39). This expression is correct up to 0{e^,5^) . 



B. the case far from the horizon 



In the region far from the horizon, spacetime will become almost fiat. In this region we assume that the rate of 
change of l/^{x) and v{x) is sufficiently small. As w e have seen for the asymptotic form of solutions in Sec. 2, we 
have four independent solutions since the ODE (3J) is of fourth order . For the solutions with short wavelength 
corresponding to fc±, we can use WKB approximation to solve Eq. (2.18). On the other hand, for the solutions with 
long wavelength corresponding to k±g, we can solve Eq. (2.18) perturbatively by treating the correction due to the 
modification of dispe rsio n relation as small. To be strict, we restrict our consideration to the region x > X2, where X2 
is that given in Eq. (4.3). In this region, we assume that the following relations 



^d I 1 d ^ 
O ' V dx 



1 



d 



(1 



.2^ < 



(1 — f 2) dx ^ 1 — ' 

are satisfied. As for higher order differentiations, we also assume that they are all restricted like 



(4.41) 



dx^ n ^ 



d^ 1 

n—- < 

dx^n ^ 
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By substituting the expansion (4.4), we find that these conditions are satisfied even in the region close to the horizon. 
Here, we define a quantity e(x) :— fcpQ(i",;2)3/2 ; which reduces to e(x) near the horizon. It will be natural to assume 

that e{x) takes its largest value in the region close to the horizon, and hence e{x) is at most of O(e^) owing to the 
restriction x > X2. In the following, we construct approximate solutions valid up to O(e^) in the sense of Sip/ip. 
We begin with considering the solutions with short wavelength. Substituting the expression 



^ _ g»/^ dx' k(x') 



(4.42) 



into Eq. (3.1), we write down the equation for k{x). Neglecting the terms on which differentiations with respect to x 
acted more than three times, 



- il - v'^)k'^ - 2vujk + oj"^ 



dx 
3fc'2 



2( ) k^ - (I - V^)k ~ VLO 

koil 




(4.43) 



is obtained, where ' is used to represent a differentiation with respect to x. Denoting the left hand side of Eq. ( 4.43| ) 
by F{k), we find that the first term on the right hand side is expressed as 

2dx \ dk ) ■ 

We denote the remaining terms on the right hand side by G(k). Following the standard prescription of the WKB 
approximation, the terms which contain differentiations with respect to x is taken to be small. Accordingly, we also 
expand k{x) in accordance with the number of differentiation as 



k{x) = /fc(°)(2^) + k'^^\x) + k^^\x) + ■■■. 
After a slightly long but a straight forward calculation, we obtain 

1 + 



(4.44) 



±kony/l~ t;2 ( 1 ± ve{x) - 



-€'ix) + Oie^) 



(1 + 2^2)^2 



-In 



1-W2 ^ 2fcon(l-w2)5/2 

(±kon{l - v^f/\l ± Ave{x) + 0(e2))) 



1 d 
2dx 

(±kon{l-v'f/' + Avw + ---), 

2 dx \ / 



= ±kon\/l - ^2. 



e^jx) 



Aiv^v'^ + (1 - "2) (u{vv'y + isvv'n 



+ (l-i;2)2 An 



n 



n 



(4.45) 



(4.46) 



(4.47) 



Now we turn to the solutions with small absolute values, i.e., k±s. In this case, we cannot use the WKB approxima- 
tion because the wavelength is not necessarily short compared with the typical scale for the background quantities to 
change. However, for the model with the standard dispersion relation, we have exact solutions for the field equation 
k±s = ±w/(l±w). We can use them as the leading order approximation, which is solutions when we neglect the terms 
related to the modification of dispersion relation. If we substitute k±s ~ ±u;/ (1 ± v) into the neglected terms, we find 
that all of them have relative order higher than e^. At fir st gl ance, one may think that the terms corresponding to 
the second and third terms in the square brackets in Eq. (4.43) give a correction of 0(6*^), but they mutually cancel 
out. As a result, the equation to determine the correction 6k±s := k±s =F t^/(l i v) is obtained as 



idx{l — v'^)6k±s ± 2uj6k±s 



2fco 



kj^l {x )clx' 



d^-d -d +d -d -S2 



s'T '^i^ F±(x). (4.48) 
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The right hand side consists of the terms related to the modification of dispersion relation, and they are small of 
Lok^l X O(e^). Different from the case for short wavelength modes, the equation to determine the correction becomes 
a differential equation. Therefore, we can say that the correction stays of O(e^) only when we are interested in the 
behavior of the solution within a small region such as xi > x > X2- Once an extended region is concerned, there is no 
reason why the correction stays of O(e^). In fact, we need to know the behavior of the solution both at infinity and 
in the matching region xi > x > X2. In such a case, the correction much larger than O(e^) can appear as explained 
in detail in Appendix 

Nevertheless, the origin of this correction is the effect of scattering due to the modified dispersion relation. Even if 
the observed spectrum of the emitted particles deviates from the thermal one due to this effect, it is still possible to 
adopt the interpretation that the spectrum is modified by the scattering during the propagation to a distant observer 
though it was initially thermal. Hence, we think that this effect should be discussed separately from the present issue. 



However, to be precise, we consider the case that the condition (4.41) replacing ^ with <C is satisfied. This is the 
case when u) is sufficiently large or when the functions v{x) and \/^l(x) rapi dly co nverge to some constants at x ^ Xj . 



In such cases, we can think of the first term on the left hand side of Eq. (4.48) as small. Then, solving Eq. (4.4S 
iteratively, we find that the correction stays of O(e^). Therefore, we obtain 

fc±.«±-^(l + 0(e2)). (4.49) 
lit; 

Consequently, we find that the solutions which behave like Q^i^i.^-^^^)^ at infinity are given by 

i'±{x) - e' r ''^'^^^''y, ij±six) = r (4.50) 

where the integral constants are chosen appropriately. Here we recall that what we wish to know is not k{x) but 
pdyk{y). Although we are keeping track of the error in the expression of k{x), we cannot evaluate the error in 
the integral k{y)dy when it is integrated from oo to the matching region where xi > x > X2- To overcome this 
difficulty, we need to make use of the existence of a conserved current 

j = Aik{x))e-^ r '"^'^'-y^'"'' , (4.51) 

where 

.,2 



A{k{x)) = LJV + {1- w )fc(^) 
~k, 



2^ {2fc(fl)(fc^i?) ~ kfj)) + Akt^R)d^ki^i) + 2k(^i^d^.k^Rf ~ dlkt^R)^ 



{h)) i-^hnAi)+0.k,R,) + ^ [l (1)') k,R„ (4.52) 

and /c(fl) and /c(/) are the real and the imaginary parts of k, respectively. The derivation of j is given in Appendix 

We evaluate this conserved current j at 2; ^ oo, where all terms that contain differentiations with respect to x 
vanish there. Adopting the normalization = 1 at a; — > 00, j is determined as 

2 

j = ujv 00 + {I - vl,)k^ - -j—kl^. (4.53) 

For I = ±,±s, by substituting kioo into the expression of j in place of k^o, we also define the conserved current ji 
corresponding to fc; . 

Owing to the conservation of j, 



By using this improved expression, we can calculate the explicit form of ^i{x) in the region xi > x > X2 without 
any ambiguity except for the constant phase factor that does not alter the absolute magnitude of the wave function. 



Expanding the expression (4.54) in powers of 5i with the substitution of Eqs. (4.4), we evaluate il^iix) keeping the 
terms up to 0{5^). Here in evaluating A(x), the terms that include becomes higher order in e, and hence we can 
neglect them all. Consequently, we obtain 
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'ip± (x) 



{2Kx)~i-^ exp ±i-koVlo\/2^xi +W± + 0{e\S^) ] , 
\ 2k 



(4.55) 



where W+ and are no different from W2 and Wi in Eq. (4.33), respectively. As noted above, there appears an 
integration constant ai that cannot be determined by the present analysis, but it is guaranteed to be a real number. 
Here, we did not give t he expl icit f orm of ip-s because we do no t use it later. 

By comparing (4.40) with (4.55), we find that the solution (4.40) obtained in the region close to the horizon is 
matched to the solutions obtained in the outer region like 



i^{x) — \/2nK e 



Vj+J 



, /"rTUJ\ ^ I ILU 



VJ+'s 



where a'_j_ are also real constants. From this expression, we can read the coefficients c+s, c_ as 



(4.56) 



c,. = -2sinh(^)r yl(l + 0(e^^^)) 

\ K / \ K / y J+s 



'27rKe 2k 



(4.57) 



The factor lo'vg appears in the fo rmu la ( 3.11 ) for the expectation value of the created particles. By differentiating 
the dispersion relation at infinity (3.4) with respect to fcoo, this factor is easily calculated as 



t^'(fcoo)Wg(fcoo) - Vkoo) 



dLu{ko 

dkry^ 



VooU! + (1 - 1^00)^0 



(4.58) 



and we find it to coincide with the conserved current. Thus, by considering the combination of a;'(fc+s)ug(fc_|_s)|c+sp 
the factor j^s in c+s cancels, and the same is also true for fc_. Finally the expectation value of the number of created 
particle is evaluated as 



N{uj) 



1 



(l + 0(e^<5^)) 



(4.59) 



V. ANALYTIC AND NUMERICAL STUDIES OF THE DEVIATION FROM HAWKING SPECTRUM 





e-1 


6° 


el 


e2 




^3/2 




2.-3/2 






^5/2 


X^ 










x" 
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^9/2 


X^ 


^3/2 





TABLE I. Table to explain the matching procedure. 

In the preceding section, to obtain the flux of the created particl es ob served in the asymptotic region, where v(x) 
is essentially constant, we propagated the near- horizon solution ( 4.40| ), which satisfies the appropriate boundary 
condition, to the infinity by matching it with the outer- region solutions ( [4.55 ), which are valid in the region distant 
from the horizon. As a result, we could determine the coefficients ci and we found that the thermal spectrum is 
reproduced up to 0(6^, 5^). 

To explain the matching procedure in more detail, here we present a table. In the construction of the near-horizon 
solution, the equation to be solved was expanded with respect to 5, and we obtained an equation which correctly 
determines the terms up to 0{b^^. These terms correspond to the first two lines in the above table. Then, expanding 
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them with respe ct to e by restricting our consideration to the region X2 < x < Xi, we obtained the expression ( 4.40| ) 
with ( 4.32| ) and ( 4.39 ), which contains the terms corresponding to the first 2x3 elements in the table. 

On the other hand, in the region distant from the horizon, we first considered an expansion of the solution with 
respect to e, and calculated the corrections up to 0{e^). Namely, the terms corresponding to the first three columns 
in the above table are obtained. As the next s tep , in the reg ion o f X2 < x < xi, we expanded this expression also 
with respect to S up to 0{S^) by substituting (iA) into Eq. ( 4.33 ), and we o btain ed the first 2x3 elements in the 
table. The expressions that we finally obtained are the outer-region solutions (4.55). 

Such two fold expansions in both schemes are simultaneously valid only in the region xi < x < xi-, where both t(x) 
and biix) are small. As mentioned above, as long as fcg is taken to be sufficiently large, this overlapping region always 
exists. Since both expressions obtained by using the above two different schemes are approximate solutions of the 
same equation, they must be identical if we t ake a n appropriate superposition of the four independent solutions. In 
fact, we foun d tha t the near-horizon solution ( 4.4C| ) can be written as a superposition of the outer-region solutions as 
given in Eq. (4.56). Now, let us look at the above table again. For each element in the table, we have assigned a power 
of X that the corresponding terms possess. As we mentioned above, we need to choose an appropriate superposition 
of the four independent outer-region solutions to achieve a successful matching. The coefficients which determine the 
weight of this superposition is nothing but c/ . Now we should note that the condition to determine the coefficients 
c; will be completely supplied by matching the x-independent elements. Once these coefficients are determined, the 
agreement of the other a;-dependent terms must be automatic from the consistency. The leading order x-independent 
elements consists of the terms of 0(e"(5°), and it is easy to see that the second lowest one consists of the terms of 
0(e^(5'^). This fact tells us that the possible modification of the coefficients c^g and c_ is at most of 0(t^8^), and 
hence the possible deviation from the thermal radiation starts only from this order. Thus we find 



1 



(l + 0(e2,53)). 



(5.1) 



Next we investigate the deviation from thermal radiation 5N/N := {N — Nthermai) / Nthermai in more detail. There 
are three different quantities of 0{d) as given in (4.1). We write them down as 



So = xv'\x=o ^ Kx, 6i^xK ^v"\x=o {Kx)bi, 62^xno{^] =: {Kx)b2, 



(5.2) 



where we introduced non-dimensional model parameters bi and 62. We also have various quantities of 0{S^) and of 
0{S^) consisted of higher derivatives of v and They are 



x^K ^v"'\x=o —'■ iKx)'^b3, x^ilo 



n 



=: {Kx)'^b4 



(5.3) 



and 



.^no ( - 



{Kx)^be, 



(5.4) 



respectively. Also, 63, 64, 65 and be are non-dimensional model parameters. One may suspect that terms including 
factors proportional to such as xk~'^v"'\x=o might appear among the correction terms of 0{e^5^). However, by 
repeating the same calculation that was given in Sec. 4 with extra higher order derivative terms, we can verify that 
such factors do not appear. From this notion, we can expect that the deviation from the thermal spectrum is given 

by 



SN 



(aooo + flooi^i + aoo2&2 + (7 other terms)) 



+ (ao3 + aisbi + 023^2) &3 + (ao4 + ai4&i + 024^2) &4 + a^b^ + aeb^ 



(5.5) 



where "a"s are some functions oi lo/k which are independent of model parameters. 

Now, we numerically confirm that the deviation actually starts from this order. The following results are obtained 
by using MATHEMATICA. As an example, let us consider a model given by 



1 

V = e 

2 



(5.6) 
(5.7) 
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For this model, we have k = 1, bi = 1 and 62 = 9/10, and the other parameters also do not vanish. For this fixed 
model, we numerically calculated the deviation from the thermal spectrum for various values of 1/A:q. The frequency 
Lu was fixed to I since our main interest is in the modes whose observed frequency at infinity becomes comparable 
with Hawking temperature {— k/2tt). The results of the numerical calculation are shown in Fig. 3 by the filed circles. 
The horizontal axis is log(l/A;o) and the vertical one is \ogdN/N. The data points arc fitted well by a lin er fu nction 
(the solid line) with its gradient, 1.99742, which perfectly agrees with the expectation represented by Eq.(5.5). 

At this point, one may notice that the deviation we obtained here is much larger than that given by Corley and 
Jacobson [B|, in which a model with 



V = i{tanh[(2Ka;)^]} 



1/2 



1, 



(5.8) 



and f2 = 1 was considered!**. The outstanding feature of their model is that 61 = 62 = ^3 = ^4 = &5 = &6 = 0. Hence, 
the terms of 0{e'^S^) in EqT(5.5) reduce to uooq k'^ /k^fll. If aooo = 0, all the terms of 0{e'^S^) in the deviation SN/N 
disappear, and it turns out to be 0{e'^S^). If so, the discrepancy between two calculations can be understood. To 
show that this is certainly the case, we repeated the numerical calculation for the same model that was discussed in 
Ref. j^. The resulting 6N/N calculated for various values of 1/fco were also plotted in Fig. 3 by the open squares. 
Again, the data points in logarithmic plot are fitted well by a linear function. But, this time, its gradient is 4.06935, 
which indicates that the deviation is actually caused by the terms of 0{e'^6^). 

Now, we can conclude that aooo = 0. Although this result might be interesting, we do not pursue this direction of 
study in this paper. Here, we would like to focus on another interesting aspect that is anticipated by the expression 
(5.5). With moderate values of model parameters, the deviation SN/N stays small for a sufficiently large kg. However, 
conversely, we can expect that the deviation from the thermal spectrum becomes large if we consider some extreme 
modificatio ns o f the special reference frame. Especially, when we consider the limiting case i n w hich fig 0, the 
expression (5.5) diverges. Although the approximation used to obtain the analytic expression ( ^.5[ ) is no longer valid 
in this limit, we can still expect that the resulting spectrum will significantly differ from the thermal one. As we 



mentioned below Eq.(2.10), there is a lower bound on Q {x). The possible smallest value of D, (x) is 1 — v , which is 
realized when v{x) = 0. Hence, we find v{x) ~ near a; = in this limiting case. Recall that v{x) was the coordinate 
velocity of the integration curves of u. Hence, vanishing v{x) means that we adopt the reference frame corresponding 
to the static observers. 

Here, we present the results of our numerical calculation, which shows that the deviation 6N/N can be large for 
some cases. Since we also want to demonstrate that the drastic change of spectrum can occur just in the consequence 
of the change of the special reference frame, we varies only the function 7(i), which was defined in E q.(|2.4| ). The 
model of v{x) is kept unchanged. As for a model of v{x), we assume the same form that is given in Eq.( |5.6| ) As for 
7(x), we adopt 
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(1/2) 



(5.9) 



with p E [0,1). p ~ corresponds to the original model associated with the freely falling observers, and p = 1 
corresponds to the case with = 0. With this choice of 7(x), the following two conditions are satisfied. One is that 
v{x) stays negative for all positive x. The other is that fi^ = 1. We calculated the deviation 6N/N for various values 
of p, and the results were shown in Fig. 4. As was expected, the deviation becomes large for small r2§. This plot raises 
an interesting speculation that TV might converge to in the fig — s- limit. Although we have not confirmed it yet, it 
is very likely that this is the case because the situation in this limit is very similar to the case in which we set a static 
mirror surrounding the event horizon of the black hole. The calculation for small fl^ were tried, but it was found 
to be out of validity of our present computation code. Anyway, we conclude that, even if K/kg is sufficiently small, 
the deviation from the thermal spectrum can be large if the combination n/koflo becomes large. It should be noted 
that the integration curves of u are not required to be suffered from infinite acceleration to achieve such a small but 
non-zero JIq. 



**T. Jacobson suggested us the existence of this discrepancy. 
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FIG. 3. The logarithmic plot of SN/N as a function of 1/fco for two different choice of v{x). The filled circles(our model) 
and the open squares(model in Ref.[5]) represent the numerical data points for respective models. Each solid line corresponds 
to a linear function which fits the data points. 




FIG. 4. A plot of 5N/N as a function of 1/no for the model given by Eq.(5.6). 

VI. CONCLUSION 

We studied the particle creation in a model which is a generalization of the Unruh's toy model. In his model, the 
field equation for a scalar field is modified by introducing a non-standard dispersion relation. To do so, we necessarily 
violates the Lorentz invariance. This radical change of theory is originally motivated by the possible existence of the 
effect due to the unknown physics at Planck scale. However, as is explained in Introduction, there is another point of 
view, on which it is also meaningful to study this model as an effective theory which takes into account the interaction 
between various fields even if we believe that the Lorentz invariance is exact. 

In the original model, the dispersion relation is modified on the basis of the freely falling observers. In our present 
work, we generalized the choice of the reference frame with respect to which we set the non-standard dispersion 
relation. Extending the analytic method developed by Corley Q, we have shown that the thermal spectrum of 
radiation from a black hole is almost reproduced as long as the modification of the special reference frame is not too 
extreme. In this analysis, we assumed lu fa k, where to is the frequency of the emitted photon observed at the spatial 
infinity and k is the surface gravity of the black hole. We have also obtained a strong suggesting that the deviation 
from the exact thermal spectrum appears from 0{k'^ /kg), where fcg is the typical wave number corresponding to the 
modification of the dispersion relation. This speculation has been confirmed numerically. 

Of course, we should not stress this small deviation from the Hawking spectrum. In the ordinary model with the 
Lorentz invariant dispersion relation, the thermal radiation at the temperature T for the static observer is observed as 

the thermal radiation at the temperature ^J^^T for the observer moving with the radial velocity /3. This argument 

holds in general whatever the source of the outward pointing radiation is because it is a direct consequence of = fc. 
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However, in the present modified model, the Lorentz invariance is violated from the beginning. We can easily see that 
(w — k)/u! is also of O(w^/feo). Hence, even if the exact Hawking spectrum is reproduced for one specific free-falling 
observer, it cannot be so for the other free-falling observers. 

On the other hand, the result that we obtained analytically also suggests that the deviation from the thermal 
spectrum can be large if we consider some extreme situations. With the aid of numerical methods, we also examined 
one of such extreme situations. We considered a sequence of different special reference frames which ranges from the 
case in which the observers associated with the special reference frame are freely falling into black hole to the case 
in which they are kept from falling into it. We found that, in the latter limiting case, the spectrum of radiation 
can significantly differ from the thermal one, even though io^ /k^ is small. It will be important to study the physical 
meaning of this result. But, since the central issue of this paper is to develop the analytic treatment of our new model, 
we have not performed detailed numerical studies yet. We will return to this issue in future publications. 
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APPENDIX A: THE CONTOUR MODIFIED FOR THE CORRECTION TERM 



In this appendix, we explain how we choose the contour of integration C in Eq. ( 4.17 ) in more detail. As shown 
in Ref. the contour C in Fig. || satisfies the condition that the wave function exponentially decreases inside the 
horizon when VF^^^ and higher order corrections are neglected. When W^^'^ are taken into account, however, the 
contour of integration needs to be modified. The correction W'^^\ contains the terms proportional to s^, while the 
terms of the highest power in the main component W''^^ is proportional to . As a result, |iy(^^| becomes larger 
than when |s| becomes large. Hence, from the validity of approximation, the contour of integration must be 

modified to be contained in the region that satisfies the condition |Ty("'| 3> IM^^^-*!. By comparing the absolute value 
of the term in w'^'^^ with that of the terms in I^'-^-', the allowed region for the contour to move is found to be 
restricted by 



S < Sr. 



min fcorioj ^o^oa 



(Al) 



Thus, we modify the contour not to run into infinity but to terminate at points contained in the region |s| < s„ 
shown in Fig. |l|. 



as 



Because of this modification of the contour of integration, the boundary terms in Eq. (4.7) no longer vanish. 
However, since '4'{sjnax) is exponentially small at both end points, we can expect that the correction due to the 
boundary terms is negligiblly small. 




FIG. 5. The integration contour satisfying the boundary condition that the wave function decays exponentially fast inside 
the event horizon when the terms of higher order in S can be neglected. In the directions indicated by hatched regions, ip{s) 
increases exponentially. Filled circles represent the saddle points s± and the dashed line represents the branch cut. 
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APPENDIX B: EVALUATION OF INTEGRATION ABOUT SADDLE POINTS 



In this appendix, we explain the details how to evaluate Eq. ( 4.2S| ). We first consider the exponent xf{s+). We 
evaluate it as an expansion around s = sq like 



+xfi{so) + f[{so){xsi) + -^fi{sQ){xsif - 



(Bl) 



The first term in the first line of the right hand side is zeroth order in S, and the second term vanishes identically. 
The other terms in the first line are quadratic or higher order in S. The terms in the second line are proportional to 
, ■ ■ ■, respectively. Thus we find 



xf{s+) = xfoiso) + xfiiso) + 0(6^). 
Using Eqs. ( [4.19D , ( |4.20| ) and ( [4.26| ), xf{s+) is found to be given by 



(B2) 




2 1 J 1 J 

3 " 10^° + 10-^^ 



25o 



l + \ log(A:orJov — 2kx) 



-(l + «-)+^(9 + i- 1-^(9 - 15z- 



<52 5 



Next, we evaluate f"{s+). Again we expand it around s = so as 

\f"{s+) = -J^{so) + ^r^'{so){xs,) + --- 

+-/r(so) + 4/"'(so)(x5i) + .--. 

X x^ 

The power indices with respect to 5 of the respective terms in the first line on the right hand side are 0, 1, 
in the second line are 1, 2, • • •. Hence, the expression for / {s+)/x which is correct up to 0{5) is given by 



1 



1 



1 



1 



/"(«+) « -/o (^o) + ^f'^'{so){xs,) + -r;{so) 



(B3) 



(B4) 
Those 

(B5) 



As for this factor, it is not necessary to find the second order correction in e. Hence, we can use the truncated 
expressions for xsq and xsi obtained by discarding the terms of 0(e) in Eqs. (4.26) and (4.27). Substituting these 
into (039), we find 



1 



/"(^+) 



"2~ 



l^\So 



3<52 



"2~ 

















K / 





(B6) 



Finally, we evaluate the second and third terms in the square brackets on the right hand side of Eq. (4.28). As 
before, we write f"'{s+)/x'^ as 



^fi'(.^o) + ^f[''\so)ixs,) + 



(B7) 



We evaluate the order of each term on the right hand side in this equation. Then, we find that the respective terms 
in the first line are of 0{e^), 0{e^S), 0{e^S^), Oie^S^), ■ ■ ■. Those in the second line are of Oie'^S),0{e^S'^),Oie^S^), ■ ■ ■. 
One may notice that the order in the first line does not change regularly. This is because the second term in the last 
line in Eq. ( 4.19D vanishes if it is differentiated more than four time. Since the leading term in / {s+)/x is O(e^), 
we can neglect the terms of 0{e^) or higher in Eq. (B7). Furthermore, we do not have to keep the terms of 0{6^) or 
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higher. Therefore, we have only to retain the terms /Q"(so)/a;^ and f["(sQ)/x'^. Substituting xso+ ~ — (— -^^^ 
into these two terms, (x/"'(s+))^/(a;/"(s+))^ is evaluated as 



2V 2 

As for a;/("''(s+)/(2;/"(s+))^, similarly we have 



(B8) 



(B9) 



The order of respective terms in the first line are 0(e'*), 0(e^5), • • •, and those in the second line are 0{e^6), 0{e^6^^ 
This time, only the term that we must keep is f['^\so)/x^. Therefore, we find 



Substituting all the above results into Eq. (4.28), finally we obtain Eq. (4.5(^). 



(BIO) 



APPENDIX C: WAVE PROPAGATION IN THE MODIFIED MODEL 



In this appendix, by solving Eq. ( 4.48 ) in a simple model, we show that ip+s, which becomes at x ^ oo, 

develops into a superposition of two modes given by ~ exp ij^ dx'k^^l{x') and by ~ exp ifdx'k^_^l{x') in the 

region of small x . He re, we assume dk^s ^ 1 as the condition that the expansion with respect to Sk^s be consistent. 
Formally, Eq. ( 4.48 ) can be integrated easily to obtain 



Sk+six) = 



1 — v'^{x) 



^ — ;-dx' , i ^r—TT-dx" 



H+{x'), 



(CI) 



where we introduced a constant xq, for definiteness, although the expression (CI) is independent of xq. The integration 
of (5fc+s becomes 



5k^s{x')dx' 



dx' 



^' dx"eC^^^'^"' 



H+{x") 



d 



1 



dx' \ 2uj 



Le~"J^o^^^^ / dx'e-'^o^^^^ H+{x') + — dx'H+{x'), 



2uj 



2uj 



(C2) 



where we used an integration by parts for the second equality. Thus ip^g (x) is expressed as 

/ / \ i fc',"' (x')da:' i I Sk-L-s{x')dx' 

yj+s[x) ~ e J ' e J~ 



r ''Tsi-')^-' j%k+s{x')dx' 



2uj 



dx'H+{x')] e»/'fe?'i(-')<i^' 



1^ rdx'e^: ^^^"iJ,(x') e^/^^^"i(^')-', 



(C3) 



where is a real constant. In the last step, Eq.(C!2) and — — ~ ^+2 were used. Let us denote the coefficient 

of exp ij^ dx'k'"^l{x') on the right hand side by As the probability for the waves to be scattered inward is 

proportional to it will be manifest that this scattering probability is not generally zero. 

As a simple example, let us consider the case that the spacetime is flat, i.e., v = but t-constant hypersurfaces 
can fluctuate randomly. We assume 71 (x) := 7(x) — 70 <C 1, where 70 is a x— independent constant. Furthermore, 
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we assume that fluctuations exist just in the interval between xq — A and xq + A. We assume that the fluctuations 
obeys the Gaussian random statistics characterized by 



M / dye'-'yi^.ixh.ix + y)), 



where we used ( ) to represent the ensemble average. Then, the Fourier transformation of 71 (x), 

1 



7iM := 



2n 



dxe"'''-fiix), 



satisfies 



By setting u = 0, in Eqs.dJ) and ( |2.10| ), we find 



v{x) = 7(2;), 



n{x) 



Using these equations, we obtain 
i f 2uj( — l—r- ^]dx' 



H,(x) 



Ho 



h{v)e 



4706 



with 



h{v) 



470 



(l-7o)(l + 7o)5 u 
uj , 2(-2 + 7o) , 2(3- 70) 



7iH + 0(7?), 



(C4) 

(C5) 
(C6) 
(C7) 

(C8) 



(l-7o)(l + 7o)^ (1+70)* (1 + 7)3 w 

(-8 + 370) (2-70) /t/N3 

"2(l + 7o)2 \u) 2(1 +70) \u 



(C9) 



The second term in the square bracket in Eq. (C8) does not depend on x and the contribution to [3 from this term 
can be neglected. 

Thus, we find that the coefficient [3 evaluated in the region cc < — A is given by 



lip i-x 



dx' 9{x' ~ [xo - A))e{{xo + A) - x')e ■'^0 H+{x') 



— — / dx' —J j dv h{v)^i{v)e^^''o 



2w 



xo-A 



r(x'-xo) 



2kg 



^ 2sin((^-t/)A) _ 
dvh{v)-1i{y) — e 



where we introduced Cj := 2lo/{1 — 7q). Then, with the aid of Eq. (C6), is evaluated as 



27rfc4 



dv\h{v)\'^ 



If A is sufficiently large, we can use the approximation 

sin ((J) - v)A)V 



Therefore, finally we obtain 



2uj^A 



2uj 



sin ((a) — h')A) 
(w-j/) 

=i 'kA5{v — Co). 

2uj 



n{v). 



(CIO) 



(Cll) 



(C12) 



l-7o^ 



n 
2lu 



l-7o^ 



(C13) 



(l-7o)« 

This expression is essentially proportional to (w/fco)^ x (wA). Since the scattering probability is also proportional to 
wA, the effect can be large for large A in principle. However, in reality this effect is suppressed because of the factor 
(uj/ko)*. If fco is taken to be a Planck scale, the factor {(jj/ko)* becomes extremely small, and then even the waves 
coming from the cosmological distance scale will not be affected significantly to induce some observable effects unless 
extraordinary oj is concerned. 
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APPENDIX D: DERIVATION OF THE CONSERVED CURRENT 



Here we derive the conserved current j given in Eq. (D3). First we note that ODE (3.1) can be derived from the 
variational principle of the action 



Sui — dxL 



with 



(Dl) 



£ -< [{-iuj + vdx)4>\ ■ [{i^ + vdx)(t)\ - idx4'){dx(, 



2fcQ 



1 „ 1 



— Ot — Ot 



This Lagrangian, £, is invariant under a global phase transformation of given by 



and 



using a trivial extension of the standard technique to derive the Noether current, we can show that 



dC 



d{d^x<P) 



dx 



dC 



d{dl4>) 



dL 



d{dl^) 



dx(p 



(D2) 
>. By 

(D3) 



becomes a conserved current which satisfies d^j = 0, although the present Lagrangian does not have the standard 
form in the sense that it contains d'^cj). Here we adopted the rule that the differentiation with respect to <^ or is 
performed as if and are independent. 

A pply ing the formula (D3) to the present case with the substitution (j) — e ^ k(x)dx _ g«/ (k{R)+ik^))dx ^ obtain 
Eq. (PI). 
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